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I.  INTRODUCTION 


The  MILLRACE*  high  explosive  test  sponsored  by  the  Defense  Nuclear 
Agency  (DNA)  involved  the  detonation  of  600  tons  of  an  ammonium  nitrate-fuel 
oil  (ANFO)  mixture  in  order  to  simulate  the  environmental  effects  of  the  air- 
blast  and  shock  which  would  result  from  exploding  nuclear  weapons;  .  Detona¬ 
tion  occurred  at  12:35:39  MDT  on  16  September  1981.  The  surface  explosion 
also  generated  a  large  cohesive  cloud  of  dust,  dirt  and  debris  which  rose 
from  the  detonation  site  to  an  altitude  of  about  4.4  km  mean  sea  level  (MSL), 
(2.9  km  above  ground  level,  AGL) ,  and  was  then  carried  away  from  the  blast 
site  by  the  prevailing  winds.  The  test  thus  represented  an  excellent 
opportunity  to  obtain  measurements  of  the  infrared  signature  of  heavily 
dirt-laden  clouds.  The  experiment  described  here  measured  the  infrared 
radiance  from  the  debris  cloud  as  a  function  of  wavelength,  angular  position 
relative  to  the  sun,  and  time  from  detonation.  These  measurements  could  then 
be  compared  with  similar  measurements  and  theoretical  calculations  for 
typical  water  clouds. 

In  the  IR  band  between  1.8  and  3.4  ^ ira ,  a  dirt  cloud  is  seen  by 
scattered  solar  radiation,  not  by  self-emission  as  would  be  the  case,  for 
example,  near  10  fi m.  The  background  signal  from  other  sources  in  the  field 
of  view  must  also  be  considered.  The  angular  relationships  are  shown  in 
Figure  1.  Note  in  particular  that  the  scattering  angle  is  defined  as  the 
angle  of  deviation  from  the  original  ray  path,  e.g.,  radiation  scattered 
directly  back  to  the  source  has  a  scattering  angle  of  180  deg.  In  general, 
scattering  functions  are  highly  peaked  in  the  forward  direction  (near  0  deg), 
a  phenomenon  familiar  to  everyone  who  has  seen  the  brilliant  outlining  of  the 
edges  of  a  small  cloud  which  has  moved  in  front  of  the  sun. 
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■Aircraft  Geometry  a  -  Solar  Zenith  Angle,  $  ■  Sensor  Zenith 
'  Azimuthal  Difference,  and  6  -  Scattering  Angle 


II.  EXPERIMENT  DESCRIPTION 


A.  INSTRUMENT  DESCRIPTION 

The  DEMON  instrument  is  a  small,  hand-guided  infrared  spectrometer 
designed  for  operation  on  a  Lear jet  aircraft.  A  somewhat  schematic  repre¬ 
sentation  of  the  infrared  spectrometer  is  shown  in  Figure  2.  For  the 
MILLRACE  test  the  instrument  employed  a  circular  variable  interference  filter 
(CVF)  to  scan  the  infrared  radiation  scattered  from  the  debris  cloud  in  the 
wavelength  interval  1.8  to  3.4  n m.  Prior  to  the  MILLRACE  Test  the  DEMON 
instrument  had  been  configured  as  a  6-detector  filter  radiometer  operating  in 
three  nonoverlapping  spectral  bands.  In  this  form  the  detectors  could  be 
optimized  for  the  expected  intensities  in  each  wavelength  band  though  at  the 
sacrifice  of  detailed  spectral  information.  For  MILLRACE  we  chose  to  con¬ 
centrate  on  spectral  information  and  hence  configured  the  instrument  as  an 
infrared  spectrometer. 

As  a  spectrometer,  the  DEMON  instrument  utilizes  a  circular  variable 
interference  filter  to  scan  the  region  between  1.77  and  3.42  fi m.  A  CVF 
consists  of  an  interference  filter  deposited  on  a  circular  substrate.  The 
wavelength  transmitted  by  the  CVF  varies  linearly  with  angle  around  the 
circumference  of  the  filter  so  that  as  the  filter  spins  in  the  optical  beam 
the  transmitted  wavelength  scans  from  one  end  of  the  passband  to  the  other 
and  chen  back  again.  The  percentage  bandwidth  of  the  filter  is  rela¬ 
tively  constant  at  all  wavelengths.  The  DEMON  CVF  has  an  intrinsic 
resolution  X  /  4  X  of  60.  The  optical  beam  traversing  the  filter  is 

somewhat  larger  than  optimum,  however,  resulting  in  a  slightly  degraded 
resolution  of  40  to  50.  The  CVF  is  spun  at  1.2  rad /sec  giving  a  complete 
double-scan  spectrum  in  5.2  sec.  Ribs  used  to  strengthen  the  filter  housing 
periodically  block  the  beam  and  thus  can  act  as  excellent  fiducial  marks  as 
well  as  providing  opportunities  to  determine  the  noise  level  and  zero  signal 


level  o£  Che  detectors.  The  ribs  are  placed  at  either  end  of  the  spectrum 
and  in  the  CO^  absorption  band  where  no  usable  signal  would  be  expected  in 
any  case. 

Two  ambient  temperature  PbS  detectors  manufactured  by  Santa  Barbara 
Research  Center  are  mounted  on  the  focal  plane.  The  detectors  are  ~ 1  mm 
square  and  have  response  curves  peaked  near  2.6  Jim.  The  detectors  are  cooled 
by  circulating  LNj  and  cold  N2  gas  through  cooling  tubes  in  the  mounting 
block;  the  temperature  is  monitored  by  thermocouples.  The  inst  ent  housing 
is  also  purged  with  the  cold  dry  nitrogen  gas  in  order  to  stal  ze  the  tem¬ 
perature  of  the  chopper  wheel,  avoid  condensing  water  on  tf  :ooled  com¬ 
ponents,  and  minimize  water  in  the  optical  path.  Though  non,  -  .y  ambient 
temperature  devices,  we  found  the  best  operating  temperature  for  this  pair  of 
PbS  detectors  to  be  -14°C.  This  temperature  optimizes  both  detectors  for 
the  largest  signal  and  thus,  because  we  are  not  detector-noise  or  background- 
limited,  for  the  highest  signal  to  noise  ratio  (SNR).  The  optical  signal  is 
modulated  at  500  Hz  by  a  chopper  wheel  inside  the  instrument  in  order  to  have 
a  reference  for  synchronous  detection.  The  amplifier  time  constant  is  5 
msec.  The  gains  of  the  two  detectors  are  set  electronically  to  differ  by  a 
factor  of  4,  resulting  in  an  overall  dynamic  range  of  400.  The  detectors 
view  slightly  different  portions  of  the  scene;  the  relative  position  of  each 
is  indicated  in  the  view  finder  so  that  the  experimenter  sighting  the 
instrument  can  guide  the  feature  of  interest  onto  either  of  the  detectors  or 
alternate  between  them. 

The  instrument  uses  a  single  f/1.5  doublet  lens.  The  6.95-cm  di¬ 
ameter  lens  is  made  of  germanium;  the  focal  length  is  10  cm.  The  effective 
diameter  of  the  lens  as  mounted  is  6.05  cm.  The  lens  has  an  antireflection 
coating  which  yields  ~90  percent  transmission  across  the  band  of  interest 
except  near  the  short-wavelength  cutoff  at  — 1.8  fJm.  The  focal  plane  of  the 
lens  was  determined  in  the  laboratory  with  a  laser  at  3  fi m. 


The  entire  instrument  is  mounted  inside  the  cabin  of  the  Learjet  in 
what  is  essentially  a  normal  laboratory  environment.  A  schematic  drawing  of 
the  cabin  installation  is  shown  in  Figure  3.  A  special-purpose  window 
mounting  with  two  ports  located  side-by-side  replaces  an  emergency  escape 
hatch  on  the  right  side  of  the  aircraft.  The  aft  port  is  used  by  the  DEMON 
instrument  itself  and  consists  of  two  8-in.  diameter  0.25-in.  thick  sapphire 
windows  mounted  in  a  double-pane  configuration.  An  8-mm  movie  camera  views 
through  the  left  (forward)  port  which  utilizes  plexiglass  windows,  also  in  a 
double-pane  mounting.  The  double-pane  mount  is  required  because  of  heavy 
frosting  frequently  encountered  when  operating  at  high  altitude. 
Additionally,  electrical  heating  tape  is  wrapped  around  the  periphery  of  the 
assembly . 

The  instrument  and  the  8-mm  movie  camera  are  slaved  together  on  a 
two-axis  gimballing  platform  so  that  both  instruments  always  look  in  pre¬ 
cisely  the  same  direction.  The  movie  camera  provides  both  a  view  finder  for 
spotting  the  spectrometer  at  specific  cloud  features  and  a  running 
visible-light  record  of  the  scene  at  1-sec  intervals.  A  block  of 
light-emitting  diode  digits  is  imaged  on  the  upper  edge  of  the  film  to 
provide  time  and  date.  A  hand-held  35-mm  camera  was  used  to  take  color 
photographs  of  the  general  scene;  these  photographs  were  also  time-tagged. 

Nitrogen  gas  for  cooling  the  focal  plane  and  purging  the  instrument 
is  obtained  from  the  boil-off  gas  from  a  25-liter  liquid  nitrogen  dewar. 
This  dewar  is  mounted  in  a  restraining  frame  just  below  the  instrument.  A 
rack  containing  electronics  is  mounted  just  forward  of  the  nitrogen  dewar.  A 
panel  in  this  rack  provided  the  capability  to  control  and  monitor  the  instru¬ 
ment.  Two  redundant  cassette  tape  recorders  record  the  output  from  the  two 
detectors,  the  time  code  generator,  the  aircraft  autopilot,  instrument 


gimbal-axis  position  indicators,  and  various  housekeeping  monitors.  In  addi¬ 
tion,  the  data  are  fed  to  a  Brush  recorder  so  that  signals  can  be  monitored 
in  real  time. 

B.  CALIBRATION 

Measurement  of  the  responsivity  of  each  detector  as  a  function  of 
wavelength  is  of  primary  importance  in  measuring  absolute  radiances. 
Calibrations  are  generally  performed  end-to-end,  i.e.,  all  elements  normally 
in  the  beam  in  flight  are  present  during  calibration.  The  exception  is  the 
aircraft  window;  window  transmission  over  the  entire  band  is  extremely  flat, 
however,  with  a  measured  value  of  76  percent. 

Calibration  of  the  instrument  consisted  of  measurements  of  wavelength 
versus  the  CVF  rotation  angle  and  of  the  f ield-of-view,  alignment,  and 
responsivity  as  a  function  of  wavelength  and  detector  temperature.  The 
primary  infrared  radiation  source  employed  for  the  calibration  was  a  Barnes 
Model  11-200  blackbody  source  which  operates  at  any  selected  temperature  in 
the  range  200  to  1000°C.  The  absolute  flux  calibration  was  made  with  an 
NBS  secondary  standard  blackbody  made  by  IR  Industries. 

The  angular  response  for  each  detector  is  approximately  square  and 
for  both  elevation  and  azimuth  is  ~0.6  deg  FWHM.  Small  asymmetries  appear 
as  a  result  of  known  aberrations  in  the  lens  system.  For  this  flight  the 
area  of  interest  is  typically  greater  than  the  f ield-of-view  so  that  such 
features  are  not  really  significant.  Generally  we  ignore  times  when  the 

character  of  the  scene  is  changing  rapidly  (as  a  result  of  scanning  over 

cloud  margins,  for  example). 

C.  DATA  REDUCTION 

The  data  reduction  process  used  to  obtain  spectral  signatures  of 
natural  and  artificial  clouds  is  outlined  in  Figure  A.  Obtaining  the 

spectral  signature  of  a  cloud  requires  more  than  just  processing  the 
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spectrometer  data  for  the  IR  spectra.  In  order  to  determine  the  scattering 
angle  it  is  necessary  to  know  the  elevation  and  azimuth  position  of  the  sun 
and  the  aircraft  and  of  the  spectrometer  relative  to  the  aircraft,  and  also 
the  height  of  the  cloud  and  its  distance.  The  aircraft  coordinates  are  ob¬ 
tained  from  the  ground  radar  track.  The  sun's  elevation  and  azimuth  angles 
are  readily  calculated  from  the  time  of  day  as  shown  in  Figure  5.  The 
elevation  and  azimuth  positions  of  the  spectrometer  are  obtained  from  the 
two-axis  gimballing  platform  and  are  recorded  with  the  IR  data  on  the  tape 
recorders.  The  cloud  height  and  width  are  estimated  from  photoanalysis  of 
the  ground-based  movie  films  and  our  own  Learjet  observations. 

The  processing  of  the  spectrometer's  IR  data  can  be  best  understood 
by  explaining  how  the  spectrometer  records  the  data  in  the  field.  As  the 
spectrometer  scans  a  target,  the  zero  signal  level  for  each  detector  is  also 
recorded.  This  internal  reference  is  provided  by  the  four  support  ribs  of 

the  filter  housing  which  briefly  block  the  detectors  during  each  scan.  The 

ribs  are  positioned  on  the  wheel  to  coincide  with  the  beginning  and  end  of 
the  filter  passband  and  in  the  atmospheric  transmission  minimum  between  2.6 
and  2.75  fX m.  The  signal  from  each  detector  is  synchronously  detected  and 
amplified  by  its  own  lock-in  amplifier.  Such  a  technique  is  required  to 
maximize  the  signal-to-noise  ratio  during  amplification  of  the  low-level 
signals.  Next,  the  signals  from  each  detector  are  FM  modulated  in  a  voltage 
to  frequency  conversion  and  redundantly  recorded.  Each  detector  has  a 

separate  frequency  channel.  Additional  data  such  as  the  focal  plane  tem¬ 
perature,  elevation  and  azimuth  of  the  spectrometer,  the  time,  and  the  CVF 
fiducial  mark  are  also  FM  encoded  on  separate  FM  channels  and  recorded  along 
with  the  spectrometer  data  on  both  tape  recorders.  The  clock  signal  also 
triggers  the  8-mm  movie  camera  at  1-sec  intervals.  The  35-mm  camera  is 

independently  operated  by  hand.  During  the  field  measurements,  the  crew  also 
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monitors  a  chart  recorder  which  displays  the  amplified  detector  signals  and 
the  angular  position  of  the  filter  wheel,  i.e.,  the  raw  IR  spectra.  The  crew 
base  their  selection  of  the  appropriate  optical  attenuation  for  the  spectro¬ 
meter  on  these  raw  IR  spectra. 

After  the  field  measurements,  the  data  on  the  cassette  tapes  and  an 
external  IRIG  time  code  are  transcribed  onto  a  computer-compatible  tape;  the 
data  channels  are  decoded  by  FM  discriminators  and  then  digitized  at  200 
samples/sec  into  computer  tapes.  The  200  samples /sec  digitizing  rate  is 
selected  to  match  the  5-msec  sampling  rate  (time  constant)  of  the  lock-in 
amplifier  and  to  obtain  the  maximum  real  information  available  from  the 
tape.  The  data  channels  are  unpacked  from  the  tape  and  a  correspondence 
between  the  CVF  angular  position  (equivalent  to  wavelength)  and  the  detector 
output  is  established  for  all  the  data  of  the  mission.  The  IR  spectra  are 
still  in  the  form  of  raw  voltage  levels.  To  convert  this  voltage  to  flux,  a 
300°C  blackbody  calibration  source  was  used  in  the  laboratory  to  derive  the 
instrument  response.  Using  Planck's  equation  to  determine  the  flux  and 
correcting  for  the  atmospheric  attenuation  of  the  laboratory  air  between  the 
blackbody  and  instrument,  the  spectral  response  for  this  spectrometer  was 
derived  as  shown  in  Figure  6.  Compensation  for  nonlinear  response  of  the 
electronics  is  achieved  by  comparing  any  nonlinear  output  voltage  to  the 
expected  linear  voltage  as  the  blackbody  flux  is  increased. 

During  the  data  processing  the  IR  data  are  filtered  digitally  by 
averaging  eight  consecutive  sampling  points.  This  averaging  reduces  the  data 
to  a  maximum  bandwidth  of  25  Hz.  This  decrease  in  data  bandwidth  does  not 
significantly  affect  the  instrument  resolving  power  as  it  still  allows  two 
samples  per  CVF  resolution  element  but  does  greatly  reduce  the  overall  noise 
background.  A  digital  plotter  produces  the  final  spectral  plots.  The  IR 
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Figure  6.  Typical  Spectral  Response  of  IR  Spectrometer  -  Detector  D3 


spectra,  focal  plane  temperature,  time,  and  elevation  and  azimuth  position  of 
the  spectrometer  are  also  printed. 

The  computer  program  incorporates  the  option  of  averaging  successive 
spectra  coherently.  Coherent  averaging  improves  the  signs 1-to-noise  ratio  of 
the  plotted  spectra  by  the  square  root  of  the  number  of  averaged  spectra. 
This  feature  is  particularly  important  for  backgrounds  in  which  intensity 
levels  are  low  and  the  observation  parameter  does  not  change  greatly  between 
successive  scans.  Up  to  42  seconds  of  data  (eight  double  spectra)  have  been 
averaged  to  obtain  some  of  the  data  presented  in  this  report. 

III.  FLIGHT  OBSERVATIONS 

A.  CLOUD  MOTION 

The  MILLRACE  explosive  charge  consisted  of  600  tons  of  ANFO  explosive 
material  packed  in  individual  50-lb  paper  sacks  which  were  then  stacked  on  a 
wooden  base  in  a  beehive  shape  9.1  m  in  diameter  by  11.4  m  in  height.  The 
detonation  of  the  ANFO  explosive  (approximately  equivalent  to  a  1  kt 
nuclear  blast)  produced  a  large,  heavily  dirt-laden  cloud  which  stabilized 
by  about  T+6m  at  an  altitude  of  ~4.4  km  MSL  (approximately  2.9  km 
above  the  desert  floor).  At  the  time  of  stabilization  the  cloud  was  approx¬ 
imately  2.5  km  in  diameter  and  extended  over  an  altitude  range  of  1  1/2  to  2 

km  (ignoring  the  stem).  The  total  mass  of  dirt  lofted  was  perhaps  4  to  10 

3 

kt,  based  on  the  similarly-sized  DICE  THROW  high  explosive  test  which  took 
place  some  8  km  to  the  northwest.  The  soil  at  the  site  is  a  mixture  of 
primarily  fine-grained  sand  and  silt  to  a  depth  of  3.4  [im  with  some 
coarser-grained  sand  at  greater  depth.  Because  it  had  rained  the  previous 
day  it  is  probable  that  the  debris  cloud  contained  a  higher  proportion  of 
water  than  would  normally  be  the  case  for  a  detonation  at  White  Sands.  The 
presence  of  water  clouds  at  the  stabilization  altitude  also  points  to 


16 


somewhat  greater  atmospheric  water  content  than  might  normally  be  expected  at 
this  high-desert  location. 


Figures  7  and  8  show  the  vertical  rise  of  the  peak  and  horizontal 
expansion  of  the  debris  cloud  as  a  function  of  time.  The  values  were 
determined  by  scaling  from  films  kindly  lent  to  us  by  John  Cockayne  of 
Science  Applications  Inc.  and  must  be  considered  preliminary.  The  films  used 
were  made  at  two  sites,  Atom  and  Granjean.  The  distance  and  azimuth  to  these 
sites  from  ground  zero  are  given  in  Table  1.  The  view  directions  are  almost 
perpendicular  to  one  another  yielding  a  good  picture  of  the  overall  cloud 
development.  It  was  necessary  to  pan  the  camera  to  keep  the  cloud  in  the 
field  of  view.  Normally  the  elevation  and  azimuth  of  the  camera  are 
indicated  by  a  coded  matrix  at  the  bottom  of  each  frame.  We  only  had 
available  to  us  16-mm  copies  of  the  35-mm  originals;  the  copy  process  had 
deleted  part  of  the  coding  matrix  with  the  result  that  we  can  only  follow  the 
altitude  development  for  a  limited  period,  viz.,  2  min  for  the  Atom  site  and 
7  min  for  the  Granjean  site.  As  can  be  seen  in  Figure  7,  the  debris  cloud 
had  essentially  reached  the  stabilization  altitude  by  T+7  min.  Figure  7 
shows  the  altitude  of  the  cloud  top  as  determined  by  visual  sighting  and  the 
aircraft  altimeter.  This  value,  4.4  km  MSL  (2.9  km  AGL) ,  is  somewhat  below 
the  maximum  determined  photographically  from  the  ground.  Because  we  believe 
the  altimeter  readings  to  be  essentially  correct,  it  is  possible  that  some 
rebounding  may  have  occurred  after  T+7  min  or  that  there  may  be  some  errors 
in  the  photographic  determinations.  (Note  that  there  have  been  first-order 
corrections  to  our  cloud  height  estimates  for  geometric  effects  arising  from 
wind-driven  motion  of  the  cloud.  The  wind  velocity  versus  altitude  readings, 
as  shown  in  Figure  9  were  part  of  the  ground  observations  made  by  Jack  Reed 
of  Sandia^.) 
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Figure  7.  Vertical  Rise  of  Debris  Cloud  Peak  versus  Time 
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Figure  8.  Horizontal  Expansion  of  the  Debris  Cloud  versus  Time 
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Table  1.  Photographic  Sites  Used  for  Cloud  Parameter  Determination 


Site 

Azimuthal  from 

Distance  (km) 

Ground  Eero  (deg) 

Atom 

37 

17.2 

Granjean 

291 

19.2 

According  to  Ref-  4,  the  debris  cloud  created  by  a  large  explosion 
will  rise  buoyantly  for  about  2  min,  after  which  turbulent  expansion  will 
continue  until  the  cloud  reaches  stabilization.  The  turbulent  expansion  may 
give  the  appearance  of  continued  buoyant  rise.  Internal  motion  after 
stabilization  is  driven  by  diffusion  and  wind  shear.  The  height-time  profile 
in  Figure  7  does  in  fact  show  a  break  in  the  slope  near  2  min.  The  cloud 
rose  at  about  12  m/s  for  the  first  2  min  (up  to  1.6  km  AGL)  and  then  at  about 
7  m/s  for  the  next  3  min  (up  to  3  km  AGL).  Nevertheless,  the  center  of  the 
expanding  cloud  rose  1.2  km  between  T+2  and  5  min  (from  1.0  to  2.2  km  AGL), 
indicating  that  the  major  part  of  the  rise  was  still  driven  by  buoyancy.  At 
stabilization  the  cloud  was  about  1.5  km  from  top  to  bottom. 

B.  FLIGHT  OPERATIONS 

The  flight  plan  was  structured  to  permit  characterizing  the  scattered 
radiation  as  a  function  of  wavelength,  scattering  angle  and  time  from  detona¬ 
tion.  We  recorded  the  IR  signatures  for  76  min  beginning  with  detonation  and 
ending  at  a  time  when  wind  shear  and  diffusion  had  produced  a  dispersed 
cloud.  The  water  clouds  present  throughout  the  flight  enabled  side  by  side 
comparison  of  natural  and  artificial  clouds  under  similar  environmental 
conditions.  The  major  details  of  our  flight  observations  and  maneuvers 
follow. 

The  Learjet  left  the  Alamogordo  airport  at  1200  MDT  in  preparation  for 
the  MILLRACE  test  then  scheduled  for  1230  MDT.  The  aircraft  flew  north  to  a 
designated  holding  area  some  20  km  north  of  the  actual  detonation  site. 


Safety  requirements  dictated  that  the  Learjet  not  be  present  in  the  area  over 
ground  zero  at  detonation  because  of  other  activity  in  the  airspace, 
primarily  a  light  plane  taking  photographs  and  4  instrumented  canisters  which 
were  parachuted  from  10.6  km  MSL  during  the  90  sec  prior  to  detonation.  The 
Learjet  was  therefore  to  remain  in  the  holding  area  about  20  km  north  until 
the  other  airplane  and  the  canisters,  which  were  being  tracked  by  radar,  were 
out  of  the  controlled  airspace. 

At  T-3  min  we  positioned  our  plane  in  a  straight  run  west  to  east  with 
the  starboard  observation  and  IR  instrumentation  ports  facing  toward  the  test 
site.  At  T=0  we  were  about  midway  in  our  traverse  and  at  the  minimum 
distance  of  18  km  from  ground  zero.  Our  motorized  35-mm  camera  and  IR 
instrumentation  recorded  the  first  minutes  of  the  development  of  the  debris 
cloud  (Figures  lOa-f).  We  experienced  a  very  mild,  short  turbulence  as  the 
detonation  shock  wave  passed;  the  pilots  also  observed  a  very  rapidly  expand¬ 
ing  hemispherical  rainbow  (prism)  effect  produced  by  the  shock-induced 
localized  change  in  the  air  density  and  hence  in  the  index  of  refraction.  For 
the  following  13  min  we  continued  our  observations  from  a  distance. 

At  T+10  min,  we  were  vectored  into  the  cloud  region  by  range  control. 
We  began  by  circling  the  debris  cloud  at  14,000'  MSL  (4.2  km)  in  a  right 
banking  turn  with  a  4.5  km  radius.  The  altitude  of  the  top  of  the  debris 
cloud  was  observed  to  be  at  ~4.4  km  MSL  (2.9  km  AGL) .  This  was  slightly 
lower  than  expected  based  on  nuclear  cloud  models  and  the  similarly-sized 
Misers  Bluff  II-2  test^  (see  earlier  discussion  on  cloud  heights).  After 
making  a  complete  circle  of  the  cloud  we  began  a  series  of  maneuvers  designed 
to  maximize  and  systematize  the  angular  coverage  while  emphasizing  ranges  of 
angles  of  particular  interest.  The  most  significant  of  these  maneuvers  are 
shown  in  Figures  lla-d.  The  aircraft  never  penetrated  the  debris  cloud,  nor, 
indeed,  came  closer  than  about  1  km,  but  only  circled  around  the  outside. 
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SCALK  IN  MILiS 
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GENERAL  ROAD  MAP  -  RANGE  AREA 


Figure  11a.  Flight  Path  of  Aircraft 


GENERAL  ROAD  MAP-  RANGE  AREA 


Figure  lib.  Flight  Path  of  Aircraft 
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GENERAL  ROAD  MAP -RANGE  AREA 


Figure  lid.  Flight  Path  of  Aircraft 


With  the  sun  to  the  rear,  we  first  banked  right  into  a  straight  line  heading 
of  18  deg  at  a  distance  of  ~1  km  from  the  cloud  edge  (Figure  11a),  the 

entire  revolution  around  the  cloud  forming  a  "D" -shaped  ground  track.  A 

second  revolution  was  executed  in  a  similar  manner  with  the  circular  path  on 
sunlight  side. 

At  T+20  minutes  we  rose  to  5.5  km  MSL  altitude  (  ~1.1  km  above  the  cloud 
top)  and  resumed  right  banking  turns  around  the  cloud  for  a  series  of  four 
different  bank  angles  (Figure  lib).  The  first  revolution  was  with  a  40  deg 
bank  in  a  2-mile  turn  radius;  the  second  with  a  35  deg  bank  and  2-mile 
radius;  the  third  with  a  25  deg  bank  and  2  1/2  mile  radius;  and  the  fourth 

with  a  15  deg  bank  and  4-mile  radius.  These  maneuvers  were  designed  to 

determine  the  1R  signature  variation  over  a  wide  range  of  observation  angles 
with  respect  to  the  sun. 

At  T+35  minutes,  we  returned  to  4.5  km  altitude,  essentially  level  with 
the  cloud  top,  and  concentrated  on  scanning  over  the  diffuse  edges  of  the 
cloud.  When  possible  we  would  scan  onto  nearby  water  clouds.  The  aircraft 
executed  a  continually  growing  spiral,  Figure  11c,  in  order  to  permit 
integrated  observations  to  be  made  over  larger  areas  of  the  cloud.  The 
spiral  eventually  reached  a  radius  of  about  8  km. 

By  the  time  this  sequence,  Figure  lid,  was  completed  at  T+66  min  the 
debris  cloud  had  become  relatively  diffuse,  particularly  at  its  periphery, 
but  was  still  visually  opaque  elsewhere.  In  order  to  permit  observations  at 
low  scattering  angles  we  flew  at  3.6  km  altitude  with  the  right  side  of  the 
plane  banked  upward  45  deg  toward  the  sun.  Scattering  angles  as  low  as 
perhaps  20  deg  were  achieved.  Because  the  edge  of  the  cloud  was  somewhat 
diffuse,  however,  the  scattering  conditions  may  not  be  directly  comparable 
with  readings  taken  much  earlier  in  the  flight.  The  data  gathering  was 
broken  off  at  T+76  min  by  which  time  the  dispersing  cloud  had  begun  to  merge 
visually  with  the  background  of  water  clouds. 
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c. 


EXPLOSION  GENERATED  WATER  CLOUD 


Ground  photography  shows  that  the  debris  cloud  had  developed  a  cap  of 
condensed  water  perhaps  as  early  as  T+3.5  min.  By  that  time  the  cloud  top 
had  reached  a  height  of  ~2.2  km  above  ground  (see  Figure  12).  Although  this 
water  cap  was  not  so  readily  identifiable  from  our  observing  platform  in  the 
air,  this  ground  observation  is  consistent  with  radiosonde  data  (Jack  Reed, 
private  communication^)  which  indicate  that  a  rise  of  about  5000  ft  (1.5 
km)  would  bring  the  temperature  of  the  rising  air  mass  to  the  local  dew 
point.  At  T+9.5  min  we  observed  from  the  Lear  jet  what  appeared  to  be  a 
small  convective  water  cloud  rising  from  the  top  of  the  debris  cloud.  The 
latter  had  essentially  stabilized  by  this  time.  The  small  cloud  appeared  to 
be  completely  water,  not  a  water-capped  dirt  plume,  but  it  is  not  iden¬ 
tifiable  in  the  ground  photography.  This  cloud  dissipated  after  several 
minutes.  The  cloud  (or,  more  properly,  the  plume)  is  intriguing  because  of 
the  possibility  that  it  implies  an  essentially  dirt-free  convective  column 
within  the  overall  cloud.  Green^  reported  that  a  dust-sampling  aircraft 
which  penetrated  the  debris  cloud  produced  by  DIAL  PACK  at  T+9  min  observed 
evidence  of  a  clean  air  shaft  in  the  center  of  the  cloud.  This  might  be 
related  to  the  apparent  pure  water  plume  observed  for  M1LLRACE,  even  though 
details  of  the  temporal  behavior  of  the  condensed  water  vapor  cap  differed 
considerably  for  DIAL  PACK.  (Green  attributed  the  clear  column  to  downward 
mixing  of  cool,  clean  air.)  Specifically,  DIAL  PACK  evidently  produced  a  cap 
of  condensed  water  vapor  very  early,  as  Green*’  notes  that  the  cap  started 
to  evaporate  at  about  T+3  to  5  min.  In  the  case  of  MILLRACE,  the  condensed 
water  cap  was  also  valuable  because  of  the  opportunity  it  afforded  us  to  scan 
both  water  and  dirt  plumes  at  essentially  the  same  range.  Unfortunately  that 
range  was  rather  long,  as  we  were  still  confined  to  the  holding  pattern  18  km 
north  of  the  debris  cloud  at  this  time. 


D. 


SCATTERING  GEOMETRY 


The  8un-c loud-sensor  geometry  is  illustrated  in  Figure  1.  Solar 
radiation  impinges  on  the  top  of  a  cloud  layer  at  a  zenith  angle  Ot  . 
Scattered  radiation  is  observed  by  the  sensor  which  is  positioned  on  a  ray 
with  zenith  angle  0  .  The  angle  between  the  azimuthal  directions  of  the  sun 
and  the  plane  is  defined  to  be  V.  Then  the  scattering  angle  6  is  defined  by 

cos-#  =  -(cos  y  ) (sin  Cl  sin  0  )  -  cos  Ot  cos  0 

The  probability  that  a  photon  of  light  will  be  scattered  through  a  given 
angle  9  is  generally  a  function  of  particle  size,  composition  and  shape. 

For  fixed  sun  and  sensor  zenith  angles,  Of  and  0 ,  a  limited  range  of 
scattering  angles  are  possible.  It  can  be  shown  that  minimum  0  is  given  by 
180  -|  0*  +  0  |  and  maximum  is  given  by  180-  [  a  -  |  .  For  MILLRACE  many 

of  the  cloud  observations  were  made  at  sensor  zenith  angles  near  90  deg  and 
solar  zenith  angles  around  38  deg.  Approximate  scattering  angles  were 
therefore  in  the  range  of  52  to  128  deg.  Some  observations  of  the  late 
diffuse  dust  cloud  were  made  at  angles  less  than  50  deg  and  some  observations 
throughout  ths  test  probably  reached  180  deg. 

IV.  DISCUSSION  OF  THE  OBSERVATIONS 
A.  INFRARED  OBSERVATIONS 

An  example  of  the  data  obtained  on  the  dust/debris  cloud  is  shown  in 
Figure  13.  The  plot  displays  two  successive  spectra  obtained  with  the  higher 
sensitivity  detector  during  a  single  360  deg  rotation  of  the  CVF.  Note  that 
Che  scan  begins  at  1.8/im,  moves  upward  to  3.4//m,  and  then  returns  to 
1.8  fjl m.  The  spectra  we  obtained  are  all  characterized  by  a  rather  strong 
peak  near  2.15  //m,  where  the  atmosphere  is  ~90  percent  transparent.  Strong 
atmospheric  absorptions  at  -^1.9/Um  (due  to  water  vapor)  and  from  2.4  to  2.8 
Mm  (due  to  and  CO2)  sharply  attenuate  the  flux  from  the  clouds  in 
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these  regions.  The  spectra  of  the  radiation  scattered  by  water  and  dust 
clouds  differ  most  in  the  2.8  to  3.4  fim  region,  as  can  be  seen  by  comparing 
the  spectra  in  Figures  13  and  14.  Figure  14  shows  one  of  the  spectral  pairs 
(one  rotation  of  the  CVF)  obtained  at  T+10  min  with  the  spectrometer  viewing 
on  the  water  cloud  which  emerged  above  the  dust  cloud  at  about  T+9  min  (as 
discussed  in  Section  III  C).  The  dust  cloud  can  scatter  appreciably  more 
flux  than  the  water  cloud  in  the  2.8  to  3.4- //m  region,  depending  on  particle 
size,  relative  humidity,  and  geometry. 

The  latter  variable,  geometry,  was  demonstrated  to  be  particularly 
important  during  several  circular  maneuvers  executed  during  our  flight  (see 
Section  III  A).  After  analyzing  the  data  obtained  on  three  of  these 
maneuvers,  the  following  pattern  emerged.  Each  track  began  with  a  leg  with 

an  aircraft  heading  of  ~300  deg,  or  a  scattering  angle  of  ~120  deg, 

during  which  time  little  or  no  3-  fim  signal  was  detected  from  the  dust 
cloud.  Then,  a  relatively  quick  change  to  a  more  easterly  heading 
with  0  HO  deg  was  made,  and  larger  signals  of  varying  strengths  were 

detected  at  2.8  to  3.4  /is.  During  one  such  maneuver  at  T+12  to  T+15 
minutes,  we  obtained  5  dual  spectral  scans  with  0  ~115  deg,  and  then  8  dual 

scans  with  6  ~  90  deg.  The  average  of  the  scans  for  these  two  ranges  of 

scattering  angles  are  shown  in  Figures  15  and  16.  For  this  modest  decrease 
in  scattering  angle  (4  0*  25°),  the  2.15-  fim  flux  increased  by  ~50 
percent,  while  the  3 -fim  flux  increased  by  about  a  factor  of  5.  (It  is 
possible  that  part  -  perhaps  a  factor  of  ~ 1.5  -  of  the  increase  near  3  fim 
may  be  due  to  smaller  absorption  by  water  vapor  along  the  line  of  sight  to 
the  cloud.  The  smaller  absorption  at  ~1.9  fi m  in  Figure  16  would  tend  to 
support  this  idea,  but  is  is  unlikely  that  the  differences  in  3-  fim  signal 

can  be  accounted  for  by  this  alone.  Also,  it  would  have  to  be  more  than 

coincidence  to  have  the  water  vapor  absorption  decrease  on  three  separate 
runs  just  as  the  plane's  heading  changed  in  this  manner.) 
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Figure  14.  Spectra  of  the  Explosion-Generated  Water  Cloud 


Spectra  with  a  Weak  3  pm  Signal:  0  Scattering  =  115  degrees 


As  noted  above,  these  spectra  (Figures  13  to  16)  still  exhibit  the 
effects  of  absorption  by  the  atmosphere  along  the  path  8un-cloud-observer. 
We  have  compared  the  positions  of  the  flux  peaks  near  2.97  and  3.15  f/m  with 
peaks  in  the  measured  spectra  of  stars  obtained  in  this  wavelength  from  a 
ground-based  observatory.  The  wavelengths  of  peak  signal  agree  to 
within  ~  2/3  of  a  resolution  element  in  spite  of  the  difference  in  total 
water  vapor  absorption  and  instruments  used  for  the  two  cases.  Flight 
spectra  at  various  times  and  scattering  angles  exhibit  2.97,  3.15  and  3.4-/im 
peaks.  We  attribute  these  spectral  differences  to  the  combined  effects  of 
atmospheric  absorption  along  the  line  of  sight  to  the  cloud  and  varying 
amounts  of  flux  scattered  in  the  direction  of  our  instrument. 

Table  2  shows  the  ratio,  R,  of  scattered  flux  at  2.15  jum  to  that  at 
3.15  /J m  for  a  variety  of  observing  conditions.  As  we  did  not  detect  any 
significant  3-/im  signal  from  the  small  water  cloud  above  the  dust  cloud,  even 
though  it  was  observed  at  smaller  scattering  angles  than  the  dust  cloud,  we 
empirically  can  only  derive  a  lower  limit  on  R  for  the  water  cloud.  Thus,  we 
proceeded  to  use  a  model  for  the  scattering  of  infrared  radiation  by  water 
clouds  to  derive  the  ratio  expected  for  such  a  cloud  under  our  observing 
conditions. 

B.  IR  SCATTERING  FROM  WATER  CLOUDS 

Previous  measurements  of  IR  scattering  from  water  clouds  as  well  as 
theoretical  calculations  indicated  that  the  radiance  near  3.0/um  would  be  1 
percent  or  less  of  the  2.2-^/m  radiance.  In  the  present  instrument  the  3.0-^m 
signal  would  be  at  or  below  measurement  threshold  when  the  2.2- fjm  signal  was 
in  saturation.  In  order  to  provide  an  estimate  of  the  3.0-/vm  water  cloud 
signal  for  comparison  with  the  3.0-//m  dust  cloud  signal  it  was  decided  to 
utilize  water  cloud  scattering  codes  which  had  been  developed  earlier  by  The 
Aerospace  Corporation  Chemistry  and  Physics  Department. 
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Aerospace's  Scattered  Solar  Radiance  from  Clouds,  SSRC,  code  models 
reflected  solar  radiation  from  a  cloud  as  the  solar  reflectance  from  a  uni¬ 
form  single  layer  plane-parallel,  semi-infinite,  cumulus  cloud.  The  cloud 
may  be  composed  of  spherical  water  or  ice  particles.  The  scattering 
calculation  is  based  on  the  Mie  single-scattering  model  with  an  empirical 
correction  for  multiple  scattering.  The  Aerospace  code  LOWTA  4,  an  enhanced 
version  of  the  AFGL  LOWTRAN  4  code,  provides  transmission  factors  for  the 
sun-cloud-detector  geometry.  Transmission  is  especially  reduced  in  the 
2.5-2.9-/im  band  due  to  water  vapor  and  is  generally  sensitive  to  atmospheric 
conditions.  When  it  is  run  with  radiosonde  data  from  the  M1LLRACE  test, 
LOWTA  4  provides  input  data  to  SSRC,  from  which  scattering  signatures  unique 
to  the  test  conditions  can  be  derived. 

C.  COMPARISON  OF  THEORY  AND  OBSERVATION 

The  resulting  variation  of  radiance  with  wavelength  and  scattering 
angle  (.9)  is  shown  for  a  water  cloud  in  Figures  17a  to  f.  Inband  radiance 
from  1.9  to  3.4/im  was  calculated  using  the  code.  It  is  seen  that  the 
radiance  peaks  at  smaller  0,  in  agreement  with  our  results  on  the  dust 
cloud.  However,  the  model  predicts  about  the  same  effect  at  2  fim  (a  factor 
of  1.8)  compared  to  3  [i m  (1.2)  for  scattering  by  a  water  cloud  in  direct 
contrast  with  the  results  presented  here  for  a  dust  cloud  where  the  3-  fJ  m 
signal  increases  significantly  faster  at  smaller  0  than  the  2.15-pm  signal. 
The  variation  of  radiance  with  scattering  angle  measured  at  2.7  fj m  by 
previous  experiments  (data  from  Ref.  2)  with  water  clouds  is  shown  in  Figure 
18.  Again,  the  signal  is  strongly  peaked  about  6=0  deg. 
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SCATTERING 


Figure  17e.  Theoretical  Water  Cloud  Spectra 
km  Clouds,  D  =  Surface 


RADIANT  INTENSITY,  kW/sr 


TABLE  2.  Ratio  of  Scattered  Flux  at  2.15  fxm  to  3.15  (Xm  and  to  3.4  fxm  for 

Different  Scattering  Angles 


Time 

(minutes  after 
blast) 

12. 15/13. 15 

12. 15/13. 3-3.4 

scat. 

(deg) 

Dust  Clouds 

11.4  -  11.9 

78 

13 

117 

12.0  -  12.7 

24 

12 

109 

15.0  -  15.9 

37 

13 

158 

15.9  -  16.7 

17 

12 

86 

54.3  -  54.7 

29 

18 

68 

54.7  -  55.7 

35 

15 

56 

Water  Clouds 

9.7  -  10.0 

300 

30 

98 

58.8  -  59.1 

100 

220 

127 

59.2  -  59.7 

200 

200 

130 

Note:  We  estimate  the  error  in  these  ratios  to  be  less  than  40  percent. 

Calculated  radiances  as  shown  in  Figure  17  can  be  used  to  scale  the 
observed  2.15-  /Urn  signal  on  the  water  cloud  to  obtain  an  estimate  of  the 
expected  signal  in  the  2.8  to  3.4-  fx m  region.  This  estimate  is  below  the 
level  of  detectability  with  our  present  instrument  and  is  thus  consistent 
with  our  data  on  the  water  cloud.  We  can  now  use  the  model  value  of  R  for 
the  water  cloud  (  ~200  to  ~100  depending  on  what  A  is  chosen)  to  compare  with 
the  observed  values  of  R  on  the  dust  cloud  (  ~40  to  ~10) .  There  is  thus 
about  a  factor  of  5  to  10  difference  in  the  ratio,  R,  between  the  dust  and 


water  clouds 


V.  CONCLUSIONS 


MILLRACE  proved  to  be  an  excellent  laboratory  for  performing  IR  mea¬ 
surements.  In  contrast  to  our  earlier  measurements  of  artificial  dust  clouds 
where  cloud  persistence  and  limited  scattering  angles  proved  troublesome, 
MILLRACE  provided  a  high  density,  long  persistence  cloud  of  sufficient  size 
to  facilitate  measurements  over  a  variety  of  scattering  angles  and 
altitudes.  The  explosion-generated  water  cap  along  with  natural  water  clouds 
in  the  vicinity  provided  valuable  supplemental  data.  From  the  aircraft 
platform  measurements  could  be  obtained  at  sufficient  altitude  to  reduce 
water  and  CC>2  absorption  effects. 

Interchange  of  data  among  ground  support  and  the  experimenters  provided 
valuable  information  not  available  in  our  earlier  measurements.  The  ground 
control  radar  track  was  used  in  the  calculation  of  scattering  angles.  Radio¬ 
sonde  data  from  one  experiment  were  incorporated  in  the  atmospheric  model  for 
the  LOWTA  transmission  and  SSRC  codes.  Ground-based  photographs  of  the  dust 
cloud  were  obtained  from  another  experimenter  and  were  used  to  estimate  cloud 
height  and  width. 

Measured  water  cloud  radiances  near  2.15 /im  and  radiances  obtained  from 
the  Aerospace  water  cloud  scattering  model  are  in  excellent  agreement, 
providing  experimental  verification  of  the  model  as  well  as  confirming  our 
measurement  technique.  As  was  suggested  by  earlier  experiments,  we  measured 
a  5  to  10  times  larger  radiance  in  the  3  to  3 ,4-^/m  band  for  a  dust  cloud 
compared  to  that  of  a  water  cloud.  The  difference  in  radiance  between  the 
two  types  of  cloud  is  a  function  of  cloud  height,  observer  to  cloud  distance, 
atmospheric  moisture  content  and  the  measurement  scattering  angle. 

Experience  gained  from  our  measurements  suggests  improvements  in 
instrumentation,  notably  the  choice  of  detectors  and  electronic  design,  that 


47 


would  increase  system  performance  -  especially  in  the  2.9  to  3.5-^/m  band. 
The  increased  sensitivity  would  allow  a  determination  of  water  and  dust 
cloud  radiances  over  a  wider  range  of  experimental  conditions  in  future  tests. 
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LAlOiATOtr  OfEMTIOMS 


The  Laboratory  Operations  of  The  Aerospace  Corporation  la  conducting  caper- 
laantal  and  thaoratlcal  Investigations  nacoasary  for  tha  evaluation  and  applica¬ 
tion  of  aclcntlflc  advancaa  to  nan  Military  tapaca  systems.  Versatility  and 
flexibility  have  been  developed  to  a  high  degree  by  the  laboratory  personnel  In 
dealing  with  the  many  problem*  encountered  in  the  nation's  tMpidly  developing 
apace  systens.  expertise  la  the  latest  scientific  developments  Is  vital  to  the 
accompli absent  of  tasks  related  to  these  probleaa.  The  laboratories  that  con¬ 
tribute  to  this  research  are: 

i 

Aero physics  Laboratory:  Launch  vehicle  and  reentry  aerodynamics  and  heat 
transfer ,  propulsion  cheats try  end  fluid  neehanlcs,  structural  mechanics,  flight 
dynamics;  high-teaperature  theraoaechanics,  gas  kinetics  and  radiation;  research 
in  environeental  chemistry  and  contamination;  cv  end  pulsed  chemical  laser 
development  including  chealcal  kinetics,  spectroscopy,  optical  resonators  and 
beam  pointing,  atmospheric  propagation,  laser  effects  and  counte measures. 

Chemistry  and  Physics  Laboratory:  Atnospherlc  chemical  reactions,  atno- 
spherlc  optics,  light  scattering,  state-spec If ic  chemical  reactions  and  radia¬ 
tion  transport  In  rocket  plumes,  applied  laser  spectroscopy.  laser  chealstry, 
battery  electrochemistry,  space  vacuum  and  radiation  effects  on  materials,  lu¬ 
brication  and  surface  phenoaena,  thermionic  emission,  photosensitive  aaterialt 
and  detectors,  atomic  frequency  standards,  and  bloenvlronmcntal  research  and 
monitoring. 

electronics  Research  Laboratory:  Microelectronics,  GaAs  low-noise  tr.d 
power  devices,  semiconductor  lasers,  electromagnetic  and  optical  propagation 
phenomena,  quar.tuc  electronics,  laser  communications,  Ildar,  and  electro-pptlcs; 
communication  sciences,  applied  electronics,  semiconductor  crystal  and  device 
physics,  radiometric  Imaging;  millimeter-wave  and  microwave  technology. 

Information  Sciences  Research  Officer  Program  verification,  program  trans¬ 
lation!  performance-sensitive  system  design,  distributed  architectures  fer 
spacebornc  computers,  fault-tolsrant  computer  systems,  artificial  Intelligence, 
and  microelectronics  applications. 

Materials  Sciences  Laboratory:  Development  of  new  materials:  metal  matrix 
composites,  polymers,  and  new  Forms  of  carbon;  component  failure  sns lysis  and 
reliability;  fracture  mechanics  and  stress  corrosion;  evaluation  of  materials  in 
apace  environment;  materials  performance  in  space  transportation  svsteas;  anal¬ 
ysis  of  systems  vulnarablUty  and  survivability  In  enemy-induced  envlronsents. 

tpace  Sciences  Laboratory:  Atmospheric  and  ionospheric  physics,  radiation 
from  the  atmosphere,  density  and  composition  of  the  upper  atmosphere,  aurorae 
and  alrglov;  aagnecospherlc  physics,  cosalc  rays,  generation  end  propagation  of 
plasma  waves  la  the  nagnetosphert;  solar  physics.  Infrared  astroncav;  the 
affects  of  nuclear  explosions,  magnetic  storms,-  and  solar  activity  on  the 
earth's  atmosphere,  ionosphere,  and  magnetosphere;  the  effects  of  optical, 
electromagnetic,  and  particulate  radiations  in  space  on  space  systems. 


